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ABSTRACT 



^ ■ At z = 0.1055, the gamma-ray burst GRB 031203 is one of the two nearest 

GRBs known. Using observations from the Very Large Array (VLA) and Chandra 
>■ ■ X-ray Observatory we derive sub-arcsecond locaHzations of the radio and X-ray 

afterglow of this GRB. We present near-infrared observations of the supernova SN 



^ ' 20031w, which exploded in the host galaxy of the GRB 031203. Our deep, high 

Q ■ resolution Magellan/PANIC data establish that this SN is spatially coincident 

^ ■ with the radio and X-ray localizations of the afterglow of GRB 031203 to sub- 

' arcsecond precision, and is thus firmly associated with the GRB. We use image 

CL|' differencing to subtract the bright emission from the host galaxy, and measure 

O ' the time evolution of the SN between ~ 5 and ~ 50 days after the GRB. The 

^ • resulting light curve has a shape which is quite different from that of the two 

a '■ SNe previously associated with GRBs, SN 1998bw and SN 2003dh. With SN 
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20031w securely associated with this burst, we confirm that all three GRBs with 
redshifts z < 0.3 were accompanied by SN explosions. 

Subject headings: supernovae: individual (SN 20031w) - gamma rays: bursts 

1. Introduction 

The emerging association between long-duration gamma-ray bursts (GRBs) and type 
Ic supernovae (SNe Ic) is perhaps the most significant breakthrough in our understanding 
of GRBs, and may also provide new insights into the physics of core-collapse SNe and the 
deaths of massive stars. The initial strong evidence for this connection came from the spatial 
and temporal coincidence between GRB 980425 and SN 1998bw (Galama et al. 1998). More 
recently, this picture was convincingly affirmed by the detection of SN features in the optical 
afterglow spectrum of GRB 030329 (Stanek et al. 2003; Hjorth et al. 2003; Matheson et 
al. 2003). Both of these milestone discoveries resulted from the study of the nearest GRBs 
yet identified, GRB 980425 at ^ = 0.0085 and GRB 030329 at z = 0.1685. Indeed, each 
relatively rare occurrence of a GRB at low redshift {z < 0.3) provides a unique opportunity 
for further study of the GRB-SN connection. 

GRB 031203^ was detected by IBIS on board the INTEGRAL spacecraft on 2003 De- 
cember 3, at 22:01:28 UT (Gotz et al. 2003). A fading X-ray afterglow was discovered 
by XMM-Newton on 2003 December 4 (Santos-Lleo et al. 2003; Campana et al. 2003; 

Rodrigucz-Pascual et al. 2003) and found to be consistent with the location of a radio 
transient (Frail 2003; Soderberg, Kulkarni, & Frail 2003). Prochaska ct al. (2003a; 2003b; 
2004) identified the host galaxy of the X-ray and radio transients, determined its redshift 
{z = 0.1055) and studied its properties in detail. Recently, several groups (Bcrsicr et al. 
2004; Taghaferri et al. 2004a; Thomsen et al. 2004; Cobb et al. 2004) reported optical 
photometric and spectroscopic observations of this GRB, which apparently reveal the sig- 
natures of an associated SN (designated SN 20031w; Taghaferri et al. 2004b). In spite of 
the low redshift of this event, its low Galactic latitude (less than 5° from the plane) and the 
resulting Galactic extinction {E[B — V] = 1.04; Schlegel, Finkbeiner, & Davis 1998), as well 
as its bright host galaxy, make the study of this event challenging. 



-'^The possible classification of this event as an X-ray flash (XRF) was debated in recent literature (e.g., 
Prochaska et al. 2004; Watson et al. 2004; Thomsen et al. 2004). However, the high-energy observations 
presented by Sazonov, Lutovinov, & Sunyaev (2004) conclusively show that this event does not fit any of 
the commonly used definitions of XRFs. 
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To overcome these difficulties, we have undertaken near infrared (NIR) observations. 
The results of this effort are reported in this Letter. Coordinated radio and X-ray obser- 
vations, which enabled us to probe the total energy output of this sub-energetic GRB, are 
reported elsewhere (Soderberg et al. 2004). 

2. Observations 

2.1. Neeir Infrcired imaging 

The location of the X-ray and radio transients associated with GRB 031203 was observed 
with Persson's Auxiliary Nasmyth Infrared Camera (PANIC; equipped with a 1024 x 1024 
HgCdTe Hawaii I infrared focal plane array), mounted on the Baade (Magellan I) 6.5m 
telescope, during 5 nights, in the J (1.25//) and Kg (2.16/x) bands. Due to their lower 
sensitivity and poor temporal coverage, we did not detect SN 20031w in our Kg-hand data. 
We therefore report in this Letter only the analysis of our J-band observations. We obtained 
18 or 36 dithered frames during four epochs, approximately 5, 7, 50, and 81 days after the 
GRB. A journal of the observations is given in Tabic 1. The integration time for each frame 
was 120 s. We obtained a common sky frame for each data set from the stacked image cube, 
and subtracted it from each frame. We then carried out fiat fielding, and shifted each frame 
to a common reference position to create the final image. 

2.2. Image subtraction and SN photometry 

In order to extract the faint SN signal from the bright host galaxy background, we have 
used the Common Point-spread-function (PSF) image subtraction method (CPM; Gal- Yam, 
Poznanski, & Maoz, in preparation)^. We find that extra flux is detected in our three early 
epochs (obtained 5.38, 7.39 and 50.21 days after the GRB) when compared with our latest 
epoch (day 81.11; see Fig. 1). We therefore proceed in our analysis by assuming that our 
latest image contains only a negligible amount of SN fight. This is supported by J-band 
observations presented by Cobb et al. (2004). 



^CPM compares a new image with a reference one by extracting an empirical point-spread-function (PSF) 
from both images, and convolving each image with the PSF of the other. While slightly degrading the final 
seeing (as the final PSF of both images is worse than the initial PSF of either image), this procedure produces 
output images with nominally identical PSFs, while introducing a minimal amount of noise. In particular, 
we have found that this algorithm often produces difference images which have lower subtraction residuals 
near the nuclei of bright galaxies compared to other popular image subtraction packages, and thus allows 
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To accurately measure the SN flux in each of the difference frames (Fig. 1, panels a-c), 
we have added an artificial point source, of known negative flux and with the native PSF, 
to our reference image. This results in the introduction of a positive artificial star to each of 
the difference images, which should have an identical PSF to any real variable point source 
(the native PSF of the reference image convolved with the PSF of the image from which it 
was subtracted). This source is visible in panels a-c of Fig. 1, in the lower left quadrant. 
Using aperture photometry (with the aperture equal to the FWHM of the artificial star), we 
measure the flux of SN 20031w with respect to this constant source, and derive the relative 
photometry of SN 20031w. We have identified 22 objects which appear in the 2MASS point 
source catalog and fall into the PANIC field of view. Of these, the brightest 6 are saturated 
or non-linear in our images. We have used the remaining 16 objects to tie our zero-point 
to the 2MASS photometry. We then used 10 secondary calibrators in the close vicinity of 
GRB 031203 to set the absolute calibration of the artificial reference star. From the scatter 
in calibrations derived using different secondary calibrators, we estimate a final absolute 
calibration error of 0.08 mag. Assuming this zero point, we calculate the J-band magnitudes 
of SN 20031w at 5.38, 7.39 and 50.21 days after GRB 031203, which are reported in Table 
2 and plotted in Fig. 3. Our calibration yields J-band magnitudes for the host galaxy and 
the SN which are consistent with those reported by Cobb et al. (2004). 

The main source of error in the derivation of the magnitudes reported here is the uneven 
background levels of the difference images, caused by effects such as non-perfect flat-fielding 
and sky subtraction (as seen in panels a-c of Fig. 1). In order to quantify this error, we 
have added 20 artificial point sources to each of the images, at random positions sampling 

both empty sky and high-background areas of the image, with the appropriate PSF. We then 
subtracted the reference image using CPM. We photometered the resulting difference frames 
and have iteratively adjusted the fiux of the artificial sources until their mean measured fiux 
from the difference image was equal to the flux we measured for SN 20031w at the same 
epoch. Finally, we adopted the measured scatter in the magnitudes of the artiflcial sources 
as the error in the SN magnitude at that epoch. We flnd that these errors are far larger 
than other sources of error, e.g., the Poisson errors in the aperture photometry of the bright 
artiflcial reference star we used in our photometric sequence, which are therefore neglected. 
To these errors we add in quadrature our absolute calibration uncertainty given above. 



better sensitivity to variable sources which are superposed on bright galactic background. 
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Table 1. J-band observations of SN 20031w 



Epoch 


UT date 


Days after 
GRB 031203 


Number of 
120s frames 


01 


2003 December 9, 07:07:39 


5.38 


18 


02 


2003 December 11, 07:20:05 


7.39 


18 


03 


2004 January 23, 03:10:01 


50.21 


36 X 


04 


2004 February 23, 00:45:25 


81.11 


36 



Data from two consecutive nights were combined. 



Table 2. J-band photometry of SN 20031w 



UT date 


Days after 


J magnitude 




GRB 031203 




2003 December 9, 07:07:39 


5.38 


21.66 ±0.60 


2003 December 11, 07:20:05 


7.39 


20.57 ±0.21 


2004 January 23, 03:10:01 


50.21 


20.60 ±0.21 



-6- 




Fig. 1. — NIR detection of SN 20031w. Panels a-c show 37.7" x37.7" sections of the difference 
images obtained by subtracting our last observed epoch (81 days after the GRB) from the 
first three epochs, obtained at 5.38, 7.39 and 50.21 days post burst. SN 20031w is at the 
center of each frame (indicated by crosshairs) and its brightness can be compared to an 
artificial reference star added to each image (marked in panel a) which has a constant fiux. 
The Grey levels in panels a-c are linearly scaled between the peak of the artificial reference 
star and Icr of the sky noise below the mode sky level. Comparing panels b and c, note the 
similar IR fiux detected at 7 and 50 days after the GRB. For orientation purposes, panel d 
shows a section from our reference (fourth epoch) PANIC image with the same scale and 
orientation. The host galaxy of GRB 031203/SN 20031w lies at the center of panel d. North 
is up. East to the left. 
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3. Results 

3.1. Spatial Coincidence with the Radio and X-ray Afterglow of GRB 031203 

The position of the radio afterglow from VLA observations (Soderberg et al. 2004) is 
(J2000) RA 08:02:30.1833(18), Dec -39:51:03.522(78), as referenced to the ICRS (uncertain- 
ties in the final digits are given in parentheses) . This position is derived from our first-epoch 
22 GHz observation, the highest resolution of our various radio observations, which had a 
beam size of 1.14x0.29 arcsec. 

The position of the X-ray afterglow is derived from our Chandra X-ray Observatory 
observations. We observed the afterglow with Chandra in a single 21.6 ksec exposure begin- 
ning at 21:35 UT on 22 January 2004 (mean epoch 49.1 days post-burst), with the afterglow 
position at the aimpoint of the AC1S-S3 CCD. Data were reduced in the standard manner 
using the ClAO tools (v. 3.0.2). In particular, a comprehensive "wavdetect" source detection 
analysis reveals six X-ray sources on the S3 chip with coincident objects in the GSC-2.2 cat- 
alog. Using these sources we confirm the absolute astrometry to ±0.13 arcsec precision, and 
derive a position for the afterglow of (J2000) RA 08:02:30.159, Dec -39:51:03.51, with 0.18" 
uncertainty. We then identify two X-ray sources on the S3 chip with point-like sources in 
our 9 December 2003 PANIC J-band image. By registering these two X-ray source positions 
against their J-band counterparts, we are able to locate the afterglow on the PANIC image 
to Ril.2-pixel precision. 

The resulting X-ray locafization, 0.56x0.64 arcsec in size, is shown as the red eUipse in 
Figure 2, along with the VLA radio (cyan) and PANIC NIR transient (green) locafizations, 

as transferred to the same image. All ellipses are two-sigma (95% confidence level). The 
WCS on the PANIC image is derived from the positions of 23 2MASS counterparts, and 
the uncertainty in the VLA position is dominated by the 0.095x0.121 arcsec RMS of this 
mapping. The NIR localization is derived from two of our subtracted images; the two relative 
localizations differ by (0.23, 0.46) pixels and so we adopt 0.029x0.058 arcsec as our one-sigma 
uncertainty. 

Fig. 2 shows that the location of the NIR transient is consistent with that of the afterglow 
across all wavelengths, to high precision, confirming the association that has already been 
inferred from the temporal coincidence of GRB 031203 and SN20031w. 

Moreover, this location is consistent with the peak of the host galaxy light: the NIR 
offset, which has the highest precision, is 0.15±0.23 pixels West, 0.26±0.46 pixels North. At 
z = 0.1055, this 0.04±0.06 arcsec offset corresponds to a physical distance of 0.07±0.13 kpc. 



- 8 - 




Fig. 2.— PANIC J-band image of the host galaxy of GRB 031203 from 9 December 2003, 
with locahzations of the Chandra X-ray afterglow (red), the VLA radio afterglow (cyan), 
and the PANIC NIR transient (green) superposed; the position of the central maximum of 
the host galaxy light is indicated by the white ticks, and North is up and East is to the 
left. The ellipses shown are two-sigma (95%-confidence) localizations; see text for details. 
The NIR transient ellipse is 0.46x0.92 pixels (0.06x0.11 arcsec) in size. The size of the 
VLA ellipse is 0.38x0.58 arcsec; the X-ray ellipse is 0.56x0.64 arcsec. Although transient 
light contributes to this image, its effect on the appearance of the host galaxy at this color 
stretch is negligible. The localizations of the afterglow and transient are consistent across 
all wavelengths, and locate the source within 0.2 kpc of the central maximum of the host 
galaxy light. 
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3.2. The J-band Light Curve 

Figure 3 shows the J-band hght curve of SN 20031w. Our deep observations detect the 
SN hght as early as ~ 5 days after the burst (the earhest detection of this SN), and suggest 
a steep rise in flux between the first two epochs (~ 5 and ~ 7 days after the GRB). The flux 
measurements at ~ 7 and ~ 50 days after the GRB are comparable. 

To better interpret these observations, we wish to compare our data to the J-band light 
curve of SN 1998bw, the prototype and best observed GRB-associated SN. Unfortunately, 
NIR observations of SN 1998bw are scarce - the only data available are 3 epochs of photom- 
etry and spectroscopy reported by Patat et al. (2001). We have used these NIR spectra, in 
conjunction with optical spectroscopy from the same source, to calculate the I — J colors 
of SN 1998bw, redshifted to z — 0.1055, using the methods described by Poznanski et al. 
(2002). We then construct a model for the J-band light curve of SN 1998bw as it would 
appear at 2; = 0.1055 in the following manner. Based on the observations of Galama et al. 
(1998) we construct the I-band light curve of SN 1998bw dX z = 0.1055, using the extinction 
correction given by Cobb et al. (2004). The resulting curve is shown in Fig. 3 (dashed 
curve). We then use the three epochs of J — J colors calculated above, along with the ap- 
propriate (smaller) J-band extinction correction (from Prochaska et al. 2004) to calculate 
J-band anchors for our light curve model (asterisks in Fig. 3). Finally, we assume that the 
shapes of the J and J-band light curves of SN 1998bw are similar, and adjust the J-band 
light curve upward by the mean of the three color measurements we have, yielding the solid 
curve plotted in Fig. 3. Note that using the J-band light curve shape we can obtain a decent 
fit to all three J-band anchors, and thus it appears that assuming this light curve shape is 
reasonable, at least during the declining part of the light curve, sampled by the available 
NIR photometry of SN 1998bw. 

Having exhausted the relevant observational data, we turn to models of SN 1998bw. 
Iwamoto et al. (1998) present hydrodynamical models of exploding C+0 stars which repro- 
duce the optical light curve and spectra of SN 1998bw. Synthetic spectra calculated by these 
authors extend out to 1.2/x. We use these spectra to calculate synthetic I — J colors of SN 
1998bw (as it would seem aX z — 0.1055) between ~ 10 and ~ 107 days after the GRB, as 
we have done with the real date. Errors resulting from incomplete coverage of the J-band 
by the available spectra, are calculated as in Poznanski et al. (2002). We then repeat the 
process described above, creating a new set of synthetic anchors (marked by small filled cir- 
cles in Fig. 3). The distribution of these anchor points suggests again that the I and J-band 
light curves of SN 1998bw had similar shapes. Indeed, applying the average synthetic color 
correction / — J ~ 0.9 and the appropriate extinction correction to our model /-band light 
curve of SN 1998bw, the resulting model J-band light curve of SN 1998bw (light solid curve 
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in Fig. 3) fits the synthetic anchors very well. The J band model hght curves derived from 

the sparse NIR observations (boldsolid ciu^rc) and the model spectra (light solid curve) are 
reasonably consistent with each other and support the notion that the / and J-band light 
curves of SN 1998bw had similar shape. Interestingly, the rough J — J colors of SN 20031w 
measured by Cobb et al. (2004; / — J ~ 1) appear similar to the colors we derive for SN 
1998bw aX z = 0.1055. 

Comparing our model SN 1998bw light curves with the data, an obvious discrepancy is 
immediately revealed. While our late-time point is marginally consistent with the models, 
the early-time J-band data points fall significantly below model expectations (by at least 
~ 0.9 mag, a 9(T effect). This discrepancy is enhanced if SN 20031w was intrinsically brighter 
than SN 1998bw by ~ 0.4 mag, as advocated by Thomsen et al. (2004). We also note 
that the fast early rise implied by our early data points is not reproduced by either model. 
Making SN 20031w less luminous overall would make the model inconsistent with our late 
time data point, and would conflict with the observations of Thomsen et al. (2004) and 
Cobb et al. (2004). Thus, it appears that SN 20031w had a significantly different light curve 
shape than SN 1998bw, at least in the J-band. 

4. Discussion and Conclusions 

In this Letter, we have presented J-band observations of SN 20031w. The sensitivity 
and resolution of our data enabled us to pinpoint the location of SN 20031w within its host 
galaxy, and to show it is consistent with sub-arcsecond localizations of the radio and X-ray 
afterglow of GRB 031203, thus confirming the association of these two events. The precise 
NIR localization of this event also puts it within 0.2 kpc from the host galaxy center. The 
J-band light curve of SN 20031w shows a rapid initial rise (5 — 7 days after the GRB) and 
evidence for bright emission more than 50 days after the GRB. The fast early rise of SN 
2003dh, associated with GRB 030329, has been interpreted by Woosley & Heger (2003) and 
Mazzali et al. (2003) as evidence for asymmetry in the explosion. A thorough investigation 
of this possibility will probably require an analysis of our data in combination with other 
extensive data sets of optical and NIR photometry and spectroscopy collected by other groups 
(e.g., Thomsen et al. 2004; Taghaferri et al. 2004; Bersier et al. 2004; Cobb et al. 2004). 

Cobb et al (2004) have recently reported / and J-band observations of this event, 
obtained with the SMARTS 1.3m telescope. A direct comparison between our observations 
and the SMARTS data, obtained on numerous epochs, is complicated by the fact that 
these authors do not present the light curve of SN 20031w. Instead, they plot the temporal 
evolution of the combined light of the SN and its bright host galaxy, derived from aperture 
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Fig. 3. — J-band photometry of SN 20031w. Our measurements (Stars with 2a error bars, 
see text) are compared to model hght curves of SN 1998bw. The bold dashed line shows 
the J-band light curve of SN 1998bw, redshifted to z = 0.1055, assuming an extinction of 
Ai = IA mag (Cobb et al. 2004). The only available J-band spectra of SN 1998bw (Patat 
et al. 2001) are used to calculate I — J color offsets, correcting for the different extinction 
values in the / and J bands (from Prochaska et al. 2004; see text). These points are plotted 
as asterisks, and serve as our guidance in constructing a model J-band light curve for SN 
1998bw. This is done by adopting the I-band light curve shape, normalized to fit these 
three points (bold solid curve). A second model, based on J-band points calculated using 
synthetic color corrections (small circles, see text) is shown as a light solid curve. Note that 
both J-band models are inconsistent with our early data points. 
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photometry, which shows considerable scatter. It is thus hard to say whether their data show 
the same early fast rise we detect. Our observation that SN 20031w had similar flux levels 
7 and 50 days after the GRB is consistent with the reported SMARTS data. Furthermore, 
comparing their /-band data with model light curves of SN 1998bw, these authors arrive 
at the conclusion that the light curve shape of SN 20031w does not resemble that of SN 
1998bw. This is in accord with our analysis of the J-band data (Cobb et al. 2004 do not 
attempt to compare their J-band data with a model of SN 1998bw). It appears that both 
our observations and the Cobb et al. (2004) data set suggest that SN 20031w had a hght 
curve quite unlike that of SN 1998bw; with a fast rise to maximum, which appears broader, 
and perhaps showing a secondary peak in the IR. 

The above discussion provides further evidence for the diversity of SNe associated with 
GRBs [see, e.g., Thomsen et al. (2004) and Lipkin et al. (2004) for recent reviews]. With 
SN 20031w, all three nearby GRBs are firmly associated with SNe, apparently supporting 
suggestions made, e.g., by Podsiadlowski et al (2004), that all long GRBs are accompanied 
by SNe. It may well be the case that the focus of future studies should now move from 
proving the association between SNe and GRBs to an attempt to characterize the properties 
of this population of SNe. Such studies may provide new clues about the progenitors and 
engines of GRBs, by requiring viable GRB models to be able to produce the large quan- 
tities of nickel derived from the SN observations, as well as valuable insights into possible 
explosion mechanisms of core-collapse SNe, which may involve GRB-like aspherical effects 
[e.g., Khokhlov et al. (1999), see Gal- Yam et al. (2004) for further discussion]. Low-redshift 
GRBs, expected to be localized in larger numbers by the upcoming SWIFT mission, as well 
as systematic studies of local core-collapse SNe (e.g., Berger et al. 2003; Soderberg et al. 
2004, in preparation; Stockdale et al. 2004) will probably shed more light on these intriguing 
questions. 
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